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Integrability criterion for the Egorov hydrodynamic type systems is presented. 
The general solution by the generalized hodograph method is found. Examples are 
given. A description of three orthogonal curvilinear coordinate nets is discussed 
■ from the viewpoint of reciprocal transformations. 
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1 Introduction 

The theory of integrable hydrodynamic type systems 

TV 

u\ = jyi{u)vi, z = i,2,...,N (1) 

was created by B.A. Dubrovin and S.P. Novikov (see [4]) and developed by S.P. Tsarev 
(see [25]). 

An integrability of the hydrodynamic type systems (1) was investigated by S.P. Tsarev 
in [25] for the distinct characteristic velocities f*, which are determined by the algebraic 
system 

det \vi{u)-vSi\ = 0. 

If (1) can be written via the Riemann invariants r*(u) [i = 1,2,...,N) in the diagonal 
form 

then the integrability condition (so-called the "semi-Hamiltonian" property) is given by 

V — V v-> — V 

where dk = d/dr^. Any semi-Hamiltonian hydrodynamic type system (1) possesses an 
infinite set of conservation laws 

dthiyC) = d^M) (2) 

and an infinite set of commuting fiows (i.e. the functions simultaneously are functions 
of X, t and r, then the compatibility conditions (m^)* = (w1)t must be fulfilled) 

N 

< = ^^}(u)4, ^ = 1,2,..., AT (3) 

j=i 

parameterized by A^ arbitrary functions of a single variable (see [25]). The algebraic 
system 

x5\ + t4(u) = K.(u), 2, A; = 1, 2, A^. (4) 

yields (in an implicit form) a general solution for (1). This is the generalized hodograph 
method. 
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Definition ([19], [20], [23]): The semi-Hamiltonian hydrodynamic type system (1) 
is said to be the Egorov, if (1) has the couple of conservation laws 

dta = d^h, dtb = d^c. (5) 

If (1) is the Egorov hydrodynamic type system, then each commuting flow (3) has the 
similar pair of conservation laws (5) (see details in [19], [20], [23]) 

dra = d.j:h, drh = d^f. (6) 

In this paper we describe a very important class of conservation laws - the Egorov 
hydrodynamic type systems (see [3], [13], [14], [15], [17], [19], [20], [23], [25], [26]), 
which have plenty applications in different areas of pure and applied mathematics, physics, 
biology and chemistry. 

In this paper we introduce the extended hodograph method for the Egorov hydro- 
dynamic type systems; we signiflcantly improve the generalized hodograph method (4) 
for the Egorov hydrodynamic type systems; if the Egorov hydrodynamic type systems 
possess Hamiltonian structure, then series of conservation laws and commuting flows 
can be found iteratively. 

The paper is organized in the following order. In the second section the Egorov hydro- 
dynamic type systems are considered with the aid of the differential geometry (conjugate 
curvilinear coordinate nets). The Egorov basic set of conservation laws and commuting 
flows is introduced. In the third section the generalized hodograph method is adopted 
for the Egorov hydrodynamic type systems. In the fourth section the extended hodo- 
graph method is presented. In the flfth section the integrability criterion of the Egorov 
hydrodynamic type systems is given. In the sixth section orthogonal curvilinear coordi- 
nate nets with symmetric rotation coefficients are considered. In the seventh section the 
Egorov hydrodynamic type systems possessing a local Hamiltonian structure are investi- 
gated. Corresponding associativity equations are derived. In the eighth section reciprocal 
transformations preserving local Hamiltonian structure are found. In the ninth section 
the Egorov three orthogonal curvilinear coordinate nets are considered. Transformations 
connecting the Egorov hydrodynamic type systems associated with local Hamiltonian 
structures are described. In the tenth section the Egorov hydrodynamic type systems 
associated with nonlocal Hamiltonian structure are briefly mentioned. 



2 Conjugate curvilinear coordinate nets 

The theory of conjugate curvilinear coordinate nets is developed by G. Darboux (see [2]). 
The linear system 

d,H, = /3,,H,, i^k (7) 
is integrable, if (i.e. if dj{diHk) = di{djHk), i ^ j ^ k) 

d,(3^, = (3^,(3,,, i^3^k. (8) 

The functions (3^^ are said to be the rotation coefficients of conjugate curvilinear coordinate 
nets, the functions are said to be the Lame coefficients. 
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Let us take some solution of the nonlinear PDE system (8). Suppose the general 
solution Hi of the linear PDE system (7) (parameterized by arbitrary functions of a 
single variable) is found. Let us take any two arbitrary solutions H(2)i and of (7), 
then one can construct the integrable hydrodynamic type system (see [25]) 

^ = ^r^, (9) 

where the Riemann invariants r* are implicit functions of "times" and t^, which usually 
are called as x and t, respectively. The general solution of this hydrodynamic type system 
is given by the generalized hodograph method (see [25]): 
Theorem [25]: The algebraic system 

xH(^i), + tH^2)i = H„ 'i = l,2,...,N (10) 

yields a general solution (in an implicit form) of the hydrodynamic type system (9). 

Thus, the general solution of the nonlinear PDE system (8) together with the general 
solution of the linear PDE system (7) describe all possible semi-Hamiltonian hydrody- 
namic type systems (9). 

Taking arbitrary solutions tp^^'' of the adjoint linear system (see (7)) 

= PhA, 2 = 1, 2, . . . , iV, (11) 

one can construct conservation laws written in the potential form 

dzf^ = ^a^(r)rfr (12) 

7 

for the hydrodynamic type systems 

ri, = ^rl, z = l, 2, . . . ,iV, (13) 

where 

a,a^ = (14) 

It means, for example, that the hydrodynamic type system (9) has an infinite number of 
the conservation laws 

parameterized by arbitrary functions of a single variable. 

Let us consider the hydrodynamic type system (9) together with its M — 2 nontrivial 
commuting flow (see (13)) 

rl^ = ^ri„ m = 2,3,...,N. (15) 

Then the generalized hodograph method (see (10)) yields the general solution (see [14]) 

N 

xH^i)i + tH(^2)i + YJ^^^>^)^ = 

fc=3 
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where Riemann invariants r* are functions of independent variables t^. This is 
invertible transformation r^{t^,t^, ...,t^). 

In this article we restrict a consideration of the hydrodynamic type systems (9) (see 
also (13) and (15)) on the symmetric case 

P^k = Pm, ^ 7^ k. (17) 

Corresponding conjugate curvilinear coordinate nets (8) were introduced by G. Darboux 
in 1866 and investigated by D.Th. Egorov in 1901 in his thesis (see [5]). It was G. Darboux 
(see [2]) who proposed to call them (see (8) and (17)) the Egorov curvilinear coordinate 
systems. From the point of view of integrability properties a remarkable progress was 
achieved by L. Bianchi in 1915 (see [1]). He found a Backlund transformation for this 
problem and established the permutability property as well as the superposition formula 
for it in the flat case, i.e. when conjugate curvilinear coordinate net becomes to be 
orthogonal 

diPik + dkPkt + l^niiPrnk = 0. (18) 

Namely, the Egorov orthogonal curvilinear coordinate nets were extensively investigated 
at that time. 

If rotation coefficients are symmetric, then the linear problems (7) and (11) co- 
incide, and corresponding conservation laws (see (12) and (14)) can be taken in the 
symmetric form too 

(Kip = ''Y^ap^dt'^ , (19) 

7 

where 

dittis^ = H(i3)iH(^)i. (20) 
It means that one can introduce the function f2 determined by its second derivatives 

- o - ^ 

Then the commuting flows (15) are determined by the unique function Q. Thus, the 
hydrodynamic type system (9) has the Egorov couple of conservation laws (5) (see (20)), 
where 

diUii = dittu = H{2)iH{i)i, dia22 = H'^2)i^ (21) 

Definition ([19], [20], [23]): The conservation law density an is said to be the 
potential of the Egorov metric. 

Let us take particular solutions i?(/3)i of the linear system (7). If the rotation 
coefficients (3^^. are symmetric (see (17)), then the Egorov hydrodynamic type systems 
(15) can be written together in the symmetric potential form (19). 

In the next section we improve the generalized hodograph method (see (10), (16)) for 
the Egorov hydrodynamic type systems. 
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3 The generalized hodograph method 

Let us introduce field variables such that Sja^ = (i.e. = a^i, see (19) 

and (20)). 

Then the algebraic system (16) for the Egorov hydrodynamic type systems (15) 

dtkUn = 9tiafc„(a) (22) 

can be written in the form 

N N 

^t'^dittk = dih ^ ^t^dak = dh, (23) 

k=l k=l 

where the conservation law density h{ai, 02, a^) is determined by its derivatives dih = 
HiH(^i^i (see (20)) and parameterized by arbitrary functions of a single variable. The 
fluxes Pfc(a) of this conservation law (cf. (2)) 

dtkh{a) = dtipkisi) (24) 

can be found in quadratures 

dpk{aL) = — rfafc„(a). 

dan 

Lemma: The Egorov hydrodynamic type systems (22) possess an arbitrary commuting 
flow (see (3) and (6)) 

drtti = dtih{sL), drttk = dtipk{a), k = 2,3, A^. (25) 

Proof: can be obtained from the compatibility condition dfk{drai) = dr{dikai). 
Thus, we have the generalized hodograph method adopted for the Egorov hydrody- 
namic type systems. 

Theorem: A general solution of the Egorov hydrodynamic type systems (15) is given 
by the algebraic system (in an implicit form, see (23) j 

= (26) 
oak 

Corollary: The function Q can be expressed explicitly via the field variables and 
found in quadratures in two steps (see (19) and below) 

dilk = aks{si)d— , dil = likiajd- 



oak Oak 

In the next section we introduce the extended hodograph method for integrability of 
the Egorov hydrodynamic type systems based on this symmetric representation. 
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4 The extended hodograph method 



Any two-component system (1) is integrable by the hodograph method (see, for instance, 
[24]). Let us adopt this hodograph method for the Egorov hydrodynamic type systems. 
The hydrodynamic type system (1) can be written in the conservative form (5). This 
couple of conservation laws can be written in the potential form 

dy = adx + bdt, dz = bdx + c(a, b)dt. 

Using the Legendre transform ^ = ax + bt — y, = bx + ct — z the above couple of 
equations can be written in the form 

d^ = xda + tdb, dC, = xdb + tdc{a, b). 

Since Xb = ta (see the first above equation), then (see the second above equation) the 
compatibility condition {C,a)b = {(,b)a 

{X + tCb)a = {tCa)b (27) 

can be written as the linear PDE equation of the second order 

$aa + Cb^ab = Ca<^bb- (28) 

Thus, the hodograph method is the transformation from the field variables (unknown 
functions) a{x,t) and b{x,t) of the quasilinear system (1) to the new field variables 
(independent variables) x{a,b) and t{a,b) of the linear system (27). 

Remark: Since (see (26)) x = ha and t = hb, then ^ = h and ^ =p- Thus, the above 
equation (28) coincides with the equation describing conservation law densities /i(a, b) of 
the hydrodynamic type system at = b^, bt = d^c^a, b). 

Example: the nonlinear elasticity equation (see, for instance, [18]) is determined by 
the function c(a). The above linear equation is reducible to the hypergeometric equation 
if c(a) = a", c(a) = In a, c(a) = e°. 

Let us consider the couple of three component Egorov hydrodynamic type system 

bt = dxu{a, b, c), by = dxv{a, b, c), (29) 

Ct = dr,v{a, b, c), Cy = dxw{a, b, c), 

which can be written in the potential symmetric form (22) 



(30) 



where y = t^. 

If functions a{x,t,y), b{x,t,y) and c{x,t,y) are common for both hydrodynamic type 
systems (29), then the compatibility condition dt{dya) = dy{dta) satisfies identically and 
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the compatibility condition dt{dyb) = dy{dtb) is identically coincides with the compatibil- 
ity condition dt{dyc) = dy{dtc). 

Lemma: // two hydrodynamic type systems (29) commute, then the relationship be- 
tween three coefficients u{a,b,c), v{a,b,c) and w{a,b,c) is given by 



UaVb + VaVc - UbVa a + ^bVc 
% = Oa , 



Uc Uc 

UcVb + - ^bVc -Ua a + VbVc s 

Ob = Oc , (31) 

Uc Uc 

^ UcVb + V'^ - UbVc - Ua ^ UaVb + VaVc - UbVa 

da = dc , 

Uc Uc 

where 

_ {UaVb + VgVc " UbVa)da + {Vg + VbVc)db + {UcVb + ~ ^bVc - Ua)dc ^^^^ 

Uc 

Under the Legendre transformation h = at^ + bt'^ + ct^ — P = bt^ + ut'^ + vt^ — 
q = ct^ + vt"^ + wt^ — , the above symmetric system (30) reduces to the three differentials 

dh = t^da + t'^db + t^dc, dp = t^db + t'^du + t^dv, dq = t^dc + t'^dv + t^dw. 

Since (see the first above equation and (26)) t^ = ha, t^ = hb, t^ = he, then (see the second 
and third above equations) the compatibility conditions {{pa)b = {Pb)a, {(la)b = {(lb)a, etc) 

it\a + t\a)b = (t^ + t^Ub + t\b)a, 
{t^Uc + t^Vc)a = {t'^Ua + t^Va)c, 
{t\lc + t\c)b = {t^ +t\b + t\b)c 

can be written as the linear PDE system (where we use the temporary notation $ = ^^) 

Ua^bb + Va^bc = ^aa + Ub^ab + Vb^ac, 

Ua^bc + Va^cc = Uc^ab + Vc^ac, (33) 
Uc^bb + Vc^bc = *ac + Ub^bc + Vb^cc- 

Thus, the above described transformation from the field variables (unknown functions) 
a{x, t, y), b{x, t, y) and c(x, t, y) of the quasilinear system (29) to the new field variables 
(independent variables) x(a, b, c), t{a, b, c) and y{a, b, c) of the linear system (33) is noth- 
ing but the extended hodograph method for a couple of three component hydrodynamic 
type systems. 
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The compatibility conditions ($cc)b = {^bc)c, {^cc)a = {^ac)c, {^ac)b = {^bc)a of the 
over-determined system (33) are equivalent (31). 

Remark: Since (see (26)) x = ha, t = and y = he, then (see (24)) = h, = p 
and = q (where we use temporary notation Fi = p and F2 = q). Thus, the above 
system (33) coincides with the system describing conservation law densities h{a, b, c) of 
both hydrodynamic type systems (29). 

Obviously, this construction easily can be extended on component case. 

5 Natural extra commuting flows. The integrability 
criterion 

Without lost of generality let us restrict our consideration on the first nontrivial three 
component Egorov hydrodynamic type system written in the conservative form 

at = bx, bt = dxu{a, b, c), q = dxv{a, b, c). (34) 

If this hydrodynamic type system is semi-Hamiltonian, then it must admit the natural 
commuting flow (cf. (29)) 

ay = Cx, by = d^v{a, b, c), Cy = d^w{a, b, c). (35) 

Indeed, since (34) is semi-Hamiltonian, then this hydrodynamic type system has an 
infinite series of conservation laws 

dth{a, b, c) = d^p{a, b, c) 

and commuting flows. Each of them can be written in the form (25) (see [23]) 

ar = dri:h{a, b, c), br = d.j,p{a, b, c), = dr,q{a, b, c), (36) 

where the commuting flow (35) has the conservation law 

dyh{a,b,c) = dj:q{a, b,c). 

The function h{a, b, c) satisfies the linear PDE system of the first order with variable 
coefficients. Thus, it is very difficult to find the function h{a, b, c) in general case. However, 
at least three particular solutions h{a,b,c) are given a priori (see (34)); i.e. = a, 
h[2) = b and /i(3) = c. The first choice is trivial, the second choice is given by (34). Thus, 
we can try to reconstruct an extra commuting flow, where the first conservation law is 

(ly Cx . 

Moreover, the compatibility condition dy{dta) = dt{dya) leads to the next conservation 
law 

by = dxv{a, b,c). 

Finally, the compatibility condition dy{dtb) = dt{dyb) leads to the criterion of an integra- 
bility for the three component Egorov hydrodynamic type systems (34). 
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Criterion of integrability: The Egorov hydrodynamic type system (34) has the 
commuting flow iff' the function w can he found in quadratures (32). 

The compatibility condition dy{dtc) = dt{dyc) leads to the same result. Thus, if the 
compatibility conditions are fulfilled, then the Egorov hydrodynamic type system (34) 
has the extra commuting fiow (35). 

Remark: Taking into account that b = = c = Zy = (see (22)) a new 
potential function Q can be introduced, where 

dQ = zdx + ^dt + ady. 

Then, the couple of hydrodynamic type system (34), (35) can be written in the form 

These equations separately can be considered as 2+1 quasilinear equations. However, in 
general case all of them are non-integrable (see, for instance, [9]). In some cases one 
of them can be integrable (by the method of hydrodynamic reductions; [9]). Then two 
other equations are reductions of higher order commuting fiows on a three component 
case. These three 2+1 equations are compatible if (31) are fulfilled. 

Obviously, the same integrability criterion can be derived for component hydrody- 
namic type system written in the conservative form and containing (5). 

Theorem: If the Egorov hydrodynamic type system (see (1) and (5)) is semi-Hamiltonian, 
then M — 2 nontrivial commuting flows can be found in quadratures. In such a case these 
hydrodynamic type systems can be written in the symmetric form (22). 

Proof: Consider the Egorov hydrodynamic type system written in the form 

dtai = (9^a2, dtak = 9^fcfc(a), k = 2, 3, N. (37) 

If this hydrodynamic type system is semi-Hamiltonian, then (A^ — 1) x {N — 1) symmetric 
matrix with the elements afcn(a) can be introduced, where afe2(a) = &fe(a) and all other 
elements are not determined yet. Since (37) is semi-Hamiltonian, then ai is the potential 
of the Egorov metric (see (9) and (21)), i.e. SjOi = Then (see (9)) diUk = H(^k)iH(i)i 

and dibk = H{k)iH(2)i. Then, indeed, all other components a^n can be found in quadratures 
(see (20)) 

, TT ZJ J m \ ' dmOikdmCln , m 

Definition: The commuting flows (15) written in the conservative form (22) are said 
to be the Egorov basic set. 

Remark: The above proof can be obtained without the Riemann invariants r'^. The 
compatibility conditions dtk{dpiam) = dfn.{dtkam) lead to the full set of relationships be- 
tween the coefficients akn, which are complicated in the variables (see, for instance, the 
above case (32)). 

The main reason for consideration of local Hamiltonian structures for the Egorov 
hydrodynamic type systems (34) is following: the general solution (26) is determined by 
the general solution (33) (in the three component case), which is parameterized by three 
arbitrary functions of a single variable. However, in general case this linear PDE system 
has variable coefficients. The Hamiltonian structure leads to a reduction of (33) to the 
linear ODE system, whose coefficients h^'^^ can be found recursively (see below). 
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6 Orthogonal curvilinear coordinate nets 



Local Hamiltonian structures of hydrodynamic type systems integrable by the generalized 
hodograph method are connected with the theory of orthogonal curvilinear coordinate nets 
(see [2]). The zero curvature condition (18) is a consequence of relationship of two linear 
problems (7) and (11) 

Hi = d,^, + 5^/5^,^^. (38) 

The existence of this first order transformation is equivalent (see [25]) the existence of 
local Hamiltonian structure for corresponding hydrodynamic type system (9) 

< = 9. (g'"?^) , (39) 



da'' ^ 

where the Hamiltonian is h = J h{a.)dx, dih = ipf^ H(]^i, the momentum density is P 
gika'a^/2, diP = xpf^H(^i)i, where 



and the fiat coordinates a are determined by its derivatives dia = ip^ the con- 

stant non- degenerate metric 

is given by 



9'' = (40) 



= 9,^f ^ + Y.[irm^t\ k=l,2,...,N. (41) 

If the rotation coefficients P^^. are symmetric (see (17)) then all above formulas simplify 
(see [25]). For instance, the zero curvature condition (18) reduces to 5/3^^ = 0, where 
S = T,dm is a shift operator. It means, that the rotation coefficients P^^. depend only on the 
differences of the Riemann invariants — r"^. The corresponding linear transformation 
(38) reduces to Hi = 6 Hi (where Hi and Hi are solutions of the linear system (7)); the 
linear system (41) reduces to SSl^'' = 0. It means, that the basic Lame coefficients 
Hi depend only on the differences of the Riemann invariants r" — r"^. Since the linear 
problems (7) and (11) coincide, then we are able to introduce the Lame coefficients with 
up/sub-indexes 

H^'^ = g''''H{s)i, H^k)i = gksH'f\ 

where the constant non- degenerate metric is given by (cf. (40)) 

-ik ^ ^ jj{t)ni jj(k)m ^ ^ ^ H(^i^mH(k)m- 

Let us introduce the adjoint flat coordinates = gikof^. Then the Hamiltonian hy- 
drodynamic type systems (39) can be written in the form 

dtai = (sik^^ ■ (42) 
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In this paper we restrict our consideration on two cases, where the the Egorov hy- 
drodynamic type system is written via flat coordinates (see (39)). The first case is 
determined by the condition that a (see (5)) is a flat coordinate Oi and 6 is a flat coordinate 
02- In general case a is an arbitrary conservation law density h. 



7 Local Hamiltonian structures. The complete inte- 
grability 

Without lost of generality for simplicity we restrict our consideration on three component 
Egorov hydrodynamic type systems (39). Since the symmetric constant matrix g^^ under 
a linear transformation of independent variables can be reduced to the diagonal or skew- 
diagonal case, we restrict our consideration of the first case (a and h are flat coordinates) 
on these two sub-cases. 

1. Let us consider the Egorov hydrodynamic type system (34) with the local Hamil- 
tonian structure 

r.dhi dhi dhi 

at = (^x-Q^, h = Ox-Q^, ct = Ox^, (43) 

where dhi/dc = b (see (5)). Then we are able to choose the extra commuting flow (35) 
written in the same Hamiltonian form 



ay Ox^^, Oy Ox^^, cy Ox^^, 



where dh2/dc = c (see (29)). 

Definition: The Egorov basic set of commuting flows (22) is said to be canonical if 
the first flat coordinate ai is a potential of the Egorov metric, and all other flat coordinates 
Ok are fluxes of the first conservation law dtkOi = dtia^, where the corresponding local 
Hamiltonian structure is given by (42). 

Since dh2/db = dhi/ da, both Egorov hydrodynamic type systems are determined by 
the sole function z[a,b) satisfying the famous associativity equation (see [3], [6], [8]) 



-abb 



^aab ^bbb y 

(45) 



where hi = bc+Zb, /i2 = c^/2 + Za and the momentum density P is given hj ho = ac+b'^ /2. 
Since the flux h of the potential a (see (5)) is a conservation law density for an arbitrary 
commuting flow (see (6)), then the shift operator d/dc is consistent with linear system 
(33) describing conservation law densities (i.e. dh/dc = h, where h and h are solutions 
of the linear system (33)). It is easy to understand, if to take into account that the first 
equation of an arbitrary commuting flow is written in the form (cf. (5) and (43)) 

Or = ^^'^ = (46) 

Theorem: The Egorov hydrodynamic type system (43) has three infinite series of 
conservation law densities, whose coefficients h^^^ {k = 0,1,2 and n = 0,1,2,...) can be 
found iteratively. 
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Proof: Taking into account (36) two other equations of an arbitrary commuting flow 
are written in the form 

K = = Cr = = d^q. 

ob oa 

Thus (see also (24)), 

dp^n^^'^ = P^n^dc + (gi'^) + zmP^n^ + Zahbh^^'^) db + [zabbP^^^ + Zaahh^n^) da, (47) 

= g^'=)rfc + {Zabbpl^^ + Zaabh^n^) db + {ZaabP^n^ + Zaaah^n^) da, 

where = a, pf^ = b, gj"^ = c; /if ^ = b, pf^ = c + Zbb, = ^ab] hf^ = c, pf^ = Zab, 
q^^ = Zaa- Then three infinite series of particular solutions by the generalized hodograph 
method (see (26)) are given in an implicit form (cf. (36)) 



x = q{a,b,c), t = p{a,b,c), y = h{a,b,c), (48) 

where 

2 oo 2 oo 2 oo 

fc=On=0 fc=On=0 k=On=0 

and akn are arbitrary constants. 

Complete integrability: In general component case the Egorov basic set (22) is 
canonical (see (39) and (42)) 

dt^an = dfi (^^^ = 5tiafc„(a). (49) 

It means, that 

The compatibility conditions dtk{dtman) = dtm{dtkan) lead to the WDVV equation (see 
[3]) 

9 7^ TT-T^ = q 



da^da^da^ da^daWa^ da^da^da^ da^da^da^ 
The canonical Egorov basic set must have infinite series of conservation laws 

= Stigg, k,s = 1,2,.. .,N, p = 0,1,2,... (50) 

and commuting flows (see [25]) 

dtsak = dti\^], k,s = 1,2,.. .,N, p = 0,1,2,... (51) 
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Since ai is a potential of the Egorov metric, then (see (6)) 



dtsai = dti 



(52) 



where h^^^ = as and t*^-, = t**, s = 1,2,.. .,N. Thus, dhf^/da^ = hfK The compatibihty 
conditions dfk{dtsai) = dts^{dtkai) lead to (see (50) and (51)) 



(p-i) 



(53) 



Thus, the conservation law densities h^s^^^ can be found in quadratures (see (52) and 
(53)) 

dh^F+^^ = q^^da\ s,k = l,2,...,N, p={),l,2,..., (54) 

where qfl = h^f\ The consistency of the conservation laws (53) with (49) lead to the 
relationships 



ip) 

s 



The substitution of (54) in r.h.s. of (55) implies the recursion relationship 



(55) 



,2l,(P+l) 



k,s 
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.dh. 



ip) 



n = 0,l,2,.... 



found by B.A. Dubrovin in [3]. Taking into account (54) the above formula leads to the 
iterative procedure (cf. (47)) 



i,k,s,m,n = 1,2, ...,N, p = 0,l,2. 



including (54). The complete integrability of the Egorov Hamiltonian hydrodynamic type 
systems (49) was proved in [25]. Then the general solution can be given in implicit form 
by the generalized hodograph method (see (26) and cf. (48)) 



da""' 



where appropriate constants. 

Example: The first choice /iq°^ = 
dynamic type system 



(0) 

a, py 



s=0 p=0 
.(0) 



b, Qq = c determines the Egorov hydro- 



d ma = d^[ac + — ] , d mh = 8,^; {be + hz^h + aZah - Zb) , d mc = c^x — + ^^ah + aZaa 
^0 \ 2 / ^0 ^0 \ 2 

Thus, first four conservation laws for the canonical Egorov basic set (43), (44) together 
with the above commuting flow can be written in the potential symmetric form 

( a b c P\ ( x\ 



d 



b u V R 
c V w S 
\P R S Q / 



d 



t 

y 



(56) 
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where all coefficients can be found in quadratures (see (47)) 

U = C + Zhh, V = Zab, W = Zaa, 



c2 

R = bc + bzbb + azab - Zb, ^ ^ ~2 ^ ^^"^ ~^ ~ ^^^^ 

Q = ac? + iP'c + a^Zaa + 'iahzab + h^Zbb — 2{aZa + bzb — z). 

2. Let us consider the Egorov hydrodynamic type system (34) with the local Hamil- 
tonian structure 

where dhi/da = b (see (5)). Then we are able to choose the extra commuting flow (35) 
written in the same Hamiltonian form 

- f) ^hl h - F) ^ - F) 

where dh2/da = c (see (29)). Since dh2/dc = dhi/db, both Egorov hydrodynamic type 
systems are determined by the sole function z{b, c) satisfying the famous associativity 
equation (see [3], [6], [8]) 

1 ~t~ Zbbc^bcc ZbbbZccc^ 

where hi = ab + Zc, h2 = ac+ Zb and the momentum density P is given by /jq = + be. 

Theorem: The hydrodynamic type system (43) is equivalent to the hydrodynamic type 
system (58) under the transformation x ^ t. 

Proof: Let us replace x t, dhi/dc a, dhi/db c, dhi/da ^ b; c ^ dhi/dd, 
a — > dhi/db, b dhi/dc. Then (58) transforms in (43). 

Remark: The transformation connecting the above associativity equations was found 
in [8]). 

3. In general case the Egorov hydrodynamic type system (39) has the couple of extra 
conservation laws (5), where a and b are not connected with the Hamiltonian structure 
(39). 

Let us consider the commuting flow (see (39) and (42)) 

dzdi = d^{dh/da'), 

where the extra conservation law is (see (5)) 

dza = dxtti. (59) 

Under the transformation (see [21]) x ^ 2;, the above commuting flow reduces to the 
Egorov hydrodynamic type system 

d,c, = d,{dh/dd), (60) 
where q = dh/da^, h = a^dh/da^ — h and the extra conservation law is (cf. (59)) 

d^{dh/dc^) = dztt, 
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where Oi = dh/dc^. Thus, in this case the potential of the Egorov metric a = Ci. At the 
same time the Egorov hydrodynamic type system (60) possesses the commuting flow 

dyCi = d^{dh/dd), 

where the extra conservation law is 

dyCi = d^C2. 

Thus, this general case is reduced to the simplest case described above. 

8 Reciprocal Transformations 

In the previous section we considered the simplest reciprocal transformation 

dy'' = a'^dt"', cr^ = const 

preserving the Egorov hydrodynamic type systems. In this section a more complicated 
reciprocal transformation is presented. Suppose we already know (see, for instance, the 
previous section) all conservation laws and commuting flows for the Egorov hydrodynamic 
type systems (22). The generalized reciprocal transformation (see [11]) contains M rows 
and M columns determined by M conservation laws and M — 1 commuting flows. The 
number M and the number (see (1)) do not correlate to each other in general case. The 
Egorov hydrodynamic type system (22) can be written in the potential symmetric form 
(19). If M ^ iV, then the generalized symmetric reciprocal transformation (cf. (19)) 

diji = aik{si)dt'' 

preserves a potential symmetric form and again transforms (22) to the Egorov hydrody- 
namic type systems. If M < A^, then the generalized reciprocal transformation 

dyi = aik{ai)dt\ i = l,2,...,M; dyi = aikdt\ i = M + l,...,N 

must contain a symmetric part including the potential a of the Egorov metric, while all 
other independent variables transform linearly (i.e. = const). 

Without lost of generality we restrict our consideration for simplicity on the reciprocal 
transformation 

dy^ = akdt\ dy^ = dt\ A; = 2, 3, A^. (61) 

Then the Egorov hydrodynamic type system (22) reduces to the similar set of the Egorov 
hydrodynamic type systems 

dyk ( — -]=dyi — , (9^fc— = 9„i ( a„fc - ] , k,n = 2,3, ...,N. 

" \ aij ai " ai \ a\ ) 

A recalculation of local Hamiltonian structures under generalized reciprocal transforma- 
tion is given in [11]. In this section we restrict our consideration on the above three 
component case (43). 
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Theorem: The local Hamiltonian structure (43) of the Egorov hydrodynamic type 
system (34) is invariant under the reciprocal transformation (61) 

dz = adx + bdt, dy = dt. (62) 

Proof: The Egorov hydrodynamic type system (34) has 4 local conservation laws 
associated with local Hamiltonian structure (43), where the momentum density P is the 
quadratic expression with respect to flat coordinates (see [23]) 

p = ac+—. (63) 

Under the above reciprocal transformation any conservation law density h reduces to 
h/a^. Thus, (63) reduces to 

P = ac+—, 

where P = —c/a, d = 1/a, c = —P/a, b = —b/a. 

Remark: The Hamiltonian density of the Egorov hydrodynamic type system (43) 

r, dh ~ dh ^ dh 

dc db dd 

is given hj h = h/a. 



9 Three orthogonal Egorov curvilinear coordinate sys- 
tems 

orthogonal Egorov curvilinear coordinate nets are described by the Bianchi-Darboux- 
Egorov-Lame system (8) 

orthogonal curvilinear coordinate nets were investigated in many publications (see, 
for instance, [1], [2], [12], [25], [27]; and plenty references therein). In this section 
we establish a new link connecting an inflnite set of orthogonal Egorov curvilinear 
coordinate systems: 

...-At''-At"-^S'-/3.';'-^!?-... (64) 

Without lost of generality we restrict our consideration for simplicity on a three component 
case only. The main advantage of the approach presented below is that all formulas are 
given via flat coordinates only (in comparison with the Riemann invariants). In such 
a case corresponding formulas can be easily used for a construction of solutions for the 
WDVV equation. 

Suppose all rotations coeflicients (3^-^ and the Lame coeflicients H(^k)i are given. It 
means that the couple of the Egorov hydrodynamic type systems (29) is given too. The 
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local Hamiltonian structure (39) can be reduced (in three component case) to (43) by a 
linear transformation of field variables [a^ — > a^a'*, cxj = const), because the symmetric 
constant non- degenerate metric g^^ is reducible to the diagonal or skew-diagonal form. 
1. Let us apply the reciprocal transformation (61) 



dT^^^ = dy 



(0) 



dy^'^ = dr^^^ 



(65) 



to the Egorov hydrodynamic type systems (29) (see also (43), (44)) written for simplicity 
in the potential symmetric form (56) 



d 



( \ / a(o) 5(0) ^(0) p(o) ^^ I ^(0) ^^ 

t(o) 



^(1) 



biO) ^(0) y(0) p(0) 

C(0) y{0) ^(0) 5.(0) 

V p(0) R(0) ^(0) g(0) J 



y 



(0) 
(0) 



where z^^^ satisfies the associativity equation (45). The hydrodynamic type systems 
(29) form the canonical Egorov basic set (43), (44) (i.e. the fluxes 6^°), c(°) of the 
first conservation law ^^{<^)a^^^ = d^(o)b^^\ dy^a^^^ = d^^c^'^^ are the corresponding fiat 
coordinates). Since a^^-* = l/a^^\ b^^^ = —b^^^/a^^\ c^^-* = —P^^^/a^^^ are fiat coordinates 
(and a^^) is a potential of the transformed Egorov metric), then the transformed canonical 
Egorov basic set (22) can be extracted from the above transformed potential symmetric 
form 



/ 3,(0) \ 



/I/ 



(0) 



(fe(0 



b(0) 



p(0) 
■ a(0) 



c(0) 
■a(0) 



„(0) &(0)p(0) d(o) b(0)c(0) 

a(0) " a{0) 2 a{0) 

p(0) 6{0)p(0) _ ^(0) {P'-°'>) _ g(0) c(0)p(0) 



,(0) 



a(0) 

c(0) 
'a(O) 



d 



M°)c(0) _ „,{0) cWpm _ 0(0) _ „,,(0) , 

a(0) ^ a(0) a(0) / 

The transformed canonical Egorov basic set must be given by (see (43), (44)) 

9j{i)c(^) = 93.{i)2^|i)^(i), C?j;(i)c(^) = (9^(i)Z^||)^(i), 



(66) 



where z^^^ satisfies the associativity equation (45). Comparing corresponding fluxes from 
the above equations 



z 



(1) 



(PW)- 



(0) 



z 



(1) 



5(0) p(0) 
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(6W)- 



(0) 



one can compute the new solution 2;^^^ of the associativity equation (45) in quadratures 

~ a(l)"o(0) ^Mi)"o(0)' 



where 



Moreover, the substitution (57) in the above differentials yields an explicit link 

of the solutions z^^^ and z^^^ of the associativity equation (45). 

Remark: The Rihaucour transformation. Transformations of local Hamiltonian struc- 
tures under generalized reciprocal transformations are described in [11]. In the simplest 
case (65) the comparison of the Lame coefficients H^^^^ and the transformed Lame coef- 
ficients -f^^ijj from (see (20)) 

leads to the Ribaucour transformation (see details in [25]) of the three orthogonal curvi- 
linear coordinate nets 

MO) o-(O) o-(O) 
_ /o(0) _ o-(l) _ 

Pjfc - /^ifc ^(0) ' ^(i)fc - ^ ^(0) • 

These rotation coefficients (5f^ and /J^^^'' are symmetric and satisfy the zero curvature 
condition (18). Thus, indeed, the Egorov three orthogonal curvilinear coordinate nets 
are preserved under the above Ribaucour transformation (see [1] and [25]). 
2. Let us rewrite the canonical Egorov basic set (66) in the form 

^tm^%aw = •^2/(1)4(1) fed)' ^xwz%^w = dy(i)C^^\ 

Thus, under the linear transformation of independent variables 

a;(2)=y(l), t(2)=tW, ^(2)^^(1) (gg) 
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the above canonical Egorov basic set reduces to the transformed canonical Egorov basic 
set 

5^(2) &(2) = 5^.(2) (c^^) + 2^5)^(2) j , 9y(2)6(2) = (9,(2)^^2)6(2)' 

(9i(2)c(2) = a,(2)z22)b(2), <9j^(2)c(2) = a,(2)2;^(])^(2), 

where 

„{2) _ ^(1) m _ (1) (2) _ (1) (2) _ (1) (2) _ , (1) (2) _ (1) 

" ~ a(i)a(i)' ~ a(l)fe(i)' ^ — ! -^6(2)6(2) — -^6(1)6(1) ' a(2)fe(2) ~ " ' ^a(2)a{2) — " • 

It means that the solution z^"^^ of the associativity equation (45) can be found in quadratures 

dz^'^^ = + z^(2'')'^2;y)j^(i), (69) 

where 

J2) _ (1) (1) ,5(1)^(1) _M) ^J2) _ (1) ^(1) +5(1)^.(1) 

Remark: The above right differential 

dG = ZbbdZab + bdzaa 

exists iff the function 2; is a solution of the associativity equation (45). This function 
G is the first such example of infinitely many expressions, which cannot be computed 
explicitly. A computation of higher conservation laws leads to similar differentials. For 
instance (see (43), (44)), 

dy{bc + Zb) = d^{cZab + G), dt + Z^ = dx[Zabic + Zbb) + hzaa - G]. 

Thus, the iterative replication of solutions for the associativity equation (45) can be 
based on the above two steps 

... ^ ^(-2) ^ ^(-1) ^ ^(0) ^ ^(1) _ .(2) ^ 

The negative direction means that the transformations (65) and (68) are applied in an 
inverse order. 

The main statement of this section: Each solution of the associativity equation (45) 
creates the canonical Egorov basic set, simultaneously, the corresponding basic Lame coef- 
ficients h\'^^ and rotation coefficients (3^^. It means, that an infinite set of solutions z^^^ 
{see (67), (69) and (70)) creates an infinite set of orthogonal curvilinear coordinate nets 
(64). Thus, a description of all solutions of the associativity equation (45) is equivalent 
to a description of three orthogonal curvilinear coordinate nets. 

This set of transformations has the same interpretation as the so-called "dressing 
method" in the soliton theory. 
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Choosing the initial solution of the associativity equation (45) z^^^ = 0, the first 
iteration yields the solution given by 

,(1) = ^ 
8a' 

The second iteration yields the same solution. The fourth iteration yields again the 
"zero" solution. However, in general case this iterative chain of transformations cannot 
be truncated. 

In this paper just a three component case was considered in details, but without any 
restrictions this approach can be applied for component case. 

Remark: The first transformation (67) is exactly the "inversion /" transformation 
for the WDVV equation found in [3] (see the formulas B.13, Appendix B). The second 
transformation is exactly the Legendre type transformation for the WDVV equation found 
in [3] (see the formulas B.l and B.2, Appendix B). 



10 Nonlocal Hamiltonian structures 

Theory of nonlocal Hamiltonian structures for hydrodynamic type systems was con- 
structed by E.V. Ferapontov in [7] (see also [16]). The simplest nonlocal Hamiltonian 
structure is associated with the metric of constant curvature (see also [10], [22]). The 
existence of the Hamiltonian structure for the Egorov hydrodynamic type systems (22) 
leads to the symmetry operator (see (46) for the local Hamiltonian structure). In this 
section without lost of generality we restrict our consideration for simplicity on the non- 
local Hamiltonian structure associated with the metric of constant curvature. The corre- 
sponding hydrodynamic type system can be written via special field variables in the 
conservative form (cf. (39)) 



ai = 



- eda!")^ + ea'h 



where g'^^ is a constant symmetric matrix and e is a constant curvature. 

Let us apply the reciprocal transformation (62) to the Egorov hydrodynamic type sys- 
tem (58). Then the quadratic relationship (connected with local Hamiltonian structure, 
see [23]) 

transforms in another type quadratic relationship (connected with the above nonlocal 
Hamiltonian structure; see [16] and [22]) 

1 _ 

- — aia2 — 0304, 

where ai = 1/a, a2 = -P/a, as = — 6/a, = —c/a. The coordinate ai is the potential of 
the Egorov metric, and the corresponding transformed Egorov hydrodynamic type system 
has the above nonlocal Hamiltonian structure. Thus, this nonlocal Hamiltonian structure 
is reducible to the local Hamiltonian structure (39). 
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